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The strong electromagnetic fields of heavy nuclei can produce a wide variety of two-photon and photonuclear 
reactions at relativistic ion colliders. We present recent results from the STAR collaboration on these 'ultra- 
peripheral' interactions, focusing on vector meson production and interferometry, and on e + e~ pair production. 
The vector meson interferometry occurs because of the symmetric initial state: nucleus 1 can emit a photon which 
scatters from nucleus 2, emerging as a vector meson, or vice-versa. The two processes are indistinguishable, and 
so interfere, even though the production points are separated enough that the produced mesons decay before their 
wave functions can overlap, so the system can be used for interesting tests of quantum mechanics. 



1. Introduction 

Electromagnetic interactions between relativis- 
tic heavy ions occur frequently, even at impact 
parameters b large enough that no hadronic inter- 
actions occur [1]. The electromagnetic fields may 
be treated as fields of almost-real photons, follow- 
ing the Weizsacker- Williams method; the photon 
flux scales as Z 2 , so the cross sections are large. 

These ultra-peripheral collisions (UPCs) can be 
divided into two classes, purely electromagnetic 
interactions (two-photon) and photonuclear inter- 
actions, where a photon from the field of one nu- 
cleus interacts with the other nucleus. We focus 
on the region b > 2Ra to exclude hadronic con- 
tamination. We will present 3 topics which have 
been studied by the Solenoidal Tracker at RHIC 
(STAR) collaboration: vector meson production, 
vector meson interferometry, and e + e~ pair pro- 
duction. 

Vector mesons are produced when a photon 
from the field of one nucleus fluctuates to a qq 
pair (virtual vector meson), which then scatters 
elastically from the other nucleus, emerging as a 
vector meson [2], as is shown in Fig. la. The elas- 
tic scattering is colorless, and can be described in 
terms of Pomeron exchange. For light mesons 
like the p°, the cross section rises slowly with the 
photon energy k. For small momentum transfers, 
the elastic scattering is coherent over the entire 
nuclear target. In the coherent region, the cross 



section depends on the size of the qq dipole fluc- 
tuation. 

Heavy mesons come from small dipoles with 
small interaction cross sections. Production is 
distributed evenly throughout the target, and 
a w A 2 . For sufficiently heavy mesons like the 
J/ip, the elastic scattering may be described in 
terms of multiple gluon exchange. The cross sec- 
tion is sensitive to gluon shadowing in nuclei [3] . 

Large dipoles (light mesons) interact with a 
large cross section with the first nucleon that they 
encounter. The interaction is on the surface of the 
target, so a rs A 4 / 3 . The exact scaling can be de- 
termined by a Glauber calculation [4,5]; for the 
p , the scaling is roug 

The large photon flux and the coherent scat- 
tering lead to a large cross section for p° pro- 
duction. For 200 GeV per nucleon Au-Au colli- 
sions at RHIC, the cross section for p° produc- 
tion is 590 mb, or about 8% of the hadronic cross 
section [4]. The STAR collaboration has already 
published results on p° production in Au-Au col- 
lisions at an energy of 130 GeV per nucleon [6]. 
The cross sections agreed with the Glauber cal- 
culations. 

The elastic scattering involves the strong force, 
and so has a very short range [7] . The p produc- 
tion is localized to within 1 fm of the two ions. 
There are two possibilities: either nucleus 1 emits 
a photon which scatters off nucleus 2, or vice 
versa. These two possibilities are indistinguish- 
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Figure 1. Schematic diagrams for (a) p° produc- 
tion and (b) e + e~ production. These processes 
may be accompanied by mutual Coulomb excita- 
tion, as in (c) and (d) respectively. The dashed 
lines in (c) and (d) show the factorization, Eq. 3. 



able, and are related by a parity tranformation. 
Vector mesons have negative parity, so the two 
amplitudes subtract, with a transverse momen- 
tum (pt) dependent phase factor to account for 
the separation. The cross section is [8] 

o" = \Ai - A 2 exp (ip T ■ b)\ 2 (1) 

where A\ and A 2 are the amplitudes for p° pro- 
duction from the two directions. At mid-rapidity 
A\ = A 2 and this simplifies to 

a = a [(1 - cos Op T • b)] (2) 

where do is the cross section without interference. 
The system acts as a 2-slit interferometer, with 
slit separation b [9]. Of course, b is unmeasur- 
able, and the observed pr spectrum is obtained 
by integrating Eq. (1) over b. The pt spectrum 
is suppressed for px < fi/{b), where (b) is the 
median impact parameter. 

The p decay distance, 7/3cr < 1 fm is far less 
than (b) w 46 fm. The p° decay before the two 



amplitudes can overlap, making this an interest- 
ing test of quantum mechanics [10]. 

Two-photon processes like e + e~, Fig. lb pro- 
duction are purely electromagnetic. To lowest or- 
der, it involves 2 photons, so the cross section 
scales as Z*. At RHIC, the cross section is ex- 
pected to be huge - 33,000 barns for gold beams. 
The coupling constant Za w 0.6 is large, so per- 
turbative descriptions of this process may not be 
adequate. The calculated Coulomb corrections to 
the lowest order diagram reduce the cross section 
by 25% [11]. However, recent all-orders calcula- 
tions find a cross section that matches the lowest 
order result [12]. The apparent discrepancy is not 
well understood [13]. 

In UPCs, the coupling constant Za (where a « 
1/137 is the electromagnetic coupling constant) is 
very large, so multiple interactions between a sin- 
gle ion pair are common. p° or e+e~ production 
may be accompanied by mutual Coulomb excita- 
tion, as in Figs, lc and Id. Each ion emits a 
photon, which excites the other nucleus, usually 
to a giant dipole resonance (GDR). The excited 
nucleus decays by emitting one or more neutrons. 
These neutrons serve as a distinctive experimen- 
tal signature of nuclear dissociation. The pres- 
ence of mutual Coulomb dissociation can be used 
to 'tag' collisions at moderate impact parameters, 
2R A < b <« 30 fm [14]. Except for the shared 
impact parameter, the reactions are independent 
[15]. For example, the cross section for p° pro- 
duction with mutual Coulomb dissociation is 

a = J d 2 bP p (b)P 2EXC (b)[l - P had (b)} (3) 

where P p (b) and P 2 EXc{b) are the probabilities 
for p° production and mutual excitation respec- 
tively. Here, Phad{b) is the probability of having 
a hadronic interaction; the last term has an effect 
similar to setting a minimum impact parameter 
bmin = %Ra- For gold ions, P p (2Ra) ~ 1% and 
P-iexc^R-a) ~ 30%, so the probability of multi- 
ple interactions is substantial. The multiple pho- 
ton exchange skews these interactions to smaller 
b. This leads to a harder photon spectrum, and 
as can be seen from Eq. (2), interference at larger 
Pt- 
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2. The STAR Detector 

The STAR detector studies heavy ion and po- 
larized proton collisions at the Rclativistic Heavy 
Ion Collider (RHIC) at Brookhaven National 
Laboratory. STAR is optimized to study cen- 
tral heavy ion collisions, which may contain thou- 
sands of particles [16]. However, it is also quite 
effective for studying UPCs where the final states 
contain 2-4 particles. STAR has collected UPC 
data in gold-gold collisions at energies of 130 GeV 
per nucleon (in 2000) and 200 GeV per nuclcon 
(in 2001). STAR has also studied p° photopro- 
duction in 200 GeV per nucleon deuteron on gold 
collisions [17]. In dA collisions, the photon usu- 
ally comes from the gold nucleus, removing the 
photon direction ambiguity present in the gold- 
gold collisions. This writeup will focus on the 
200 GeV per nucleon gold-gold collisions. 

Charged particles are reconstructed in a 4.2 
meter long, 4 meter diameter time projection 
chamber (TPC) [18]. The TPC is in a solenoidal 
magnetic field which has been operated at both 
0.25 and 0.5 T. At 0.5 T, the reconstruction effi- 
ciency for charged pions is high for tracks with 
transverse momentum pr > 100 MeV/c and 
pseudorapidity \r]\ < 1.15. The track position 
and specific energy loss, dE/dx were measured 
at 45 points for high momentum charged parti- 
cles with \q\ < 1. In a 0.5 T field, the dE/dx 
resolution was 8%. The TPC is surrounded by 
240 scintillator slats covering \r]\ < 1, comprising 
the central trigger barrel (CTB). Two zero degree 
calorimeters (ZDCs) are located 18 m upstream 
and downstream of the interaction point. These 
calorimeters are sensitive to neutrons from mu- 
tual Coulomb dissociation [19]. 

Two different triggers [20] were used to study 
UPCs. The 'topology' trigger selected events with 
an appropriate topology in the CTB. It divided 
the CTB into 4 quadrants: north, south, top and 
bottom, and required hits in the north and south 
quadrants. The top and bottom regions were used 
as vetoes to reject cosmic rays. 

The minimum bias trigger selected events with 
one or more neutrons in each ZDC. The ZDC sig- 
nals were required to occur within 1 nsec, restrict- 
ing these events to the central 30 cm of the TPC. 



Data from both triggers was processed identi- 
cally except that events from the CTB based 
trigger were distributed more broadly along the 
TPC axis, and consequently, were accepted in a 
broader range. 

3. p° Production 

UPC p° production has a distinctive signature - 
two oppositely charged tracks with small net pr- 
We select events with exactly 2 primary tracks 
that form a vertex. The vertexing procedure con- 
siders all of the tracks in an event, and rejects 
tracks that are inconsistent with coming from a 
single vertex. We do allow a few non-primary 
background tracks in the TPC. The topology and 
minimum bias data were treated identically ex- 
cept that allowance was made for the different dis- 
tribution of the accepted event production points. 
This study used about 1.5 million minimum bias 
and about 1.7 million topology triggers taken in a 
0.5 T field for 200 GeV per nucleon interactions. 

Figure 2a shows the pr distribution of all 
charge-0 2-track vertices. The large peak for 
Pt < 100 MeV/c is a signature of fully coher- 
ent interactions. Like-sign pion pairs (the shaded 
histogram) arc used as a background estimate; 
they are normalized to fit the unlike-sign pairs 
for pt > 250 MeV/c. This procedure treats inco- 
herent p° production (a photon scattering from a 
single nucleon in the target) as background. 

Figure 2b shows the rapidity distribution of the 
coherent p° (pairs with pr < 150 MeV/c). The 
points (data) are in excellent agreement with a 
calculation based on the soft Pomeron model and 
our detector simulation. 

Figure 2c shows the irir invariant mass, M„ 
for the sample. The data is fit to a relativis- 
tic Breit-Wigner for the p°, plus a Soding term 
to account for the interference with direct 7r + 7r~ 
production [21]. The interference shifts the peak 
of the distribution to lower M Ilr . The direct 7T7t 
to p ratio agrees with the STAR 130 GeV analysis 
[6] and with that observed by the ZEUS collabo- 
ration in 7p interactions [22]. However, because 
of the coherent enhancement, the STAR data is 
at smaller \t\ than the ZEUS analysis, so the two 
results may not be directly comparable [23]. 
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Figure 2. The (a) pr, (b) rapidity and (c) irir 
invariant mass spectra of p° using data from the 
minimum bias trigger. The mass is fit to to a 
Breit-Wigner function for the p°, plus a Soding 
term due to interference with a direct tttt com- 
ponent. The combinatoric hadronic and coherent 
e + e~ pair backgrounds are also shown. 



4. p° Interferometry 

For the interference analysis, a clean, low- 
background sample is important, so stringent 
event selection criteria were used, at some cost in 
efficiency. Events were required to have exactly 
two tracks with a vertex within 50 cm longitudi- 
nally of the center of the TPC for the minimum 
bias sample, and 100 cm for the topology sample. 
p candidates were further required to have rapid- 
ity \y\ > 0.1, to eliminate a small remaining cos- 
mic ray contamination. Finally, the pairs were re- 
quired to have 550 MeV < M„ < 920 MeV. The 
minimum mass cut removed most of the hadronic 
backgrounds, which are concentrated at low M W7r . 
These cuts reduced the like-sign pair background 
to a few percent. With these cuts, background 
from misidcntified two-photon production of lep- 
ton pairs should be very small. The sample still 
includes direct pions, which should have the same 
spin/parity and quantum mechanical behavior as 
the pion pairs from p decay. We do not distin- 
guish between the two sources. 

The interference depends on the amplitudes for 
p production on the two nuclei. Away from y = 0, 
the photon energies for the two photon directions 
differ, k% t 2 = My/2exp(±y/2); the amplitudes 
differ and the interference is less than maximal. 
Although it is not expected in the soft-Pomeron 
model, the the photon energy difference could in- 
troduce a small p° production phase difference, 
which could affect the interference [2]. This anal- 
ysis focuses on the region near mid-rapidity where 
any phase difference is small. A Monte Carlo cal- 
culation is used to find the interference for differ- 
ent rapidity ranges [4,8]. 

We use the variable t± = ip\ to study the in- 
terference. At RHIC energies, the longitudinal 
component of the 4-momentum transfer is small, 
so i ~ t±. Without interference, the spectrum 
dN/dt rs exp(— bt) for a variety of nuclear mod- 
els [8,24]. Our calculations use a Woods-Saxon 
distribution for the gold density distribution. 

Figure 3 compares the uncorrected minimum 
bias data for 0.1 < < 0.5 with two simulations, 
with and without interference. Both simulations 
include the detector response. The data has a 
significant downturn for t < 0.001 GeV 2 , consis- 
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Figure 3. Raw (uncorrected) t± spectrum for 
p° sample for 0.1 < \y\ < 0.5 for the topol- 
ogy data. The points (red) are the data. The 
dashed histogram (blue) is a simulation that as- 
sumes that there is interference, while the dot- 
dash histogram (black) is based on a calculation 
without interference. The solid black histogram 
with very few counts is the like-sign background. 



tent with the (b) = 20 fm expected for a 3-photon 
reaction [15]. This drop matches the drop seen in 
the calculation with interference, but is absent in 
the calculation without the interference. 

The efficiency corrected data are shown in Fig. 
4. Minimum bias and topology data are shown 
separately, each with two rapidity bins: 0.1 < 
\y\ < 0.5 and 0.5 < \y\ < 1.0. The efficiency is 
independent ofpx- However px smearing (resolu- 
tion) does affect the spectrum slightly. The p° px 
resolution is about 9 MeV/c, while the 1st t bin 
covers to (15 MeV/c) 2 . Interference depletes 
the first few bins, but feed down from the higher 
t bins partially repopulates them. 

The data is fit to the 3-parameter form: 

dN 

^-=aexp(-bt)[l + c(R(t)-l] (4) 
at 

where R(t) = Int(t)/Noint(t) is the ratio of the 
Monte Carlo t-spectra with and without interfer- 
ence. Here, a is the overall normalization, the 
slope b « R 2 A , and c is the degree of spectral 
modification; c = corresponds to no interfer- 



Trigger Rapidity b (GeV 2 ) c 
M.B. 0.1 - 0.5 301 ± 14 1.01 ±0.08 
M.B. 0.5 - 1.0 304 ± 15 0.78 ±0.13 
Topo. 0.1 - 0.5 361 ±10 0.71 ±0.16 
Topo. 0.5 - 1.0 368 ±12 1.22 ±0.21 
Table 1 

Results of the fits to the minimum bias (M.B.) 
and topology (topo.) data for two rapidity bins. 
The fits have x 2 of 50, 73, 81 and 50 respectively 
for 47 degrees of freedom. 



ence while c = 1 is the expected interference. 
This functional form separates the interference 
(c) from the nuclear form factor (b). 

Table 1 gives the results of the fits. At small- 
rapidities the amplitudes from the two directions 
are similar and the interference reduces the cross 
section at px = by more than at larger rapidi- 
ties. In the minimum bias data, the interference 
extends to higher px than the topology data be- 
cause the former has a smaller (b). 

The 4 c values are consistent within errors; the 
weighted average is c = 0.93 ± 0.06. The b values 
for the minimum bias and exclusive p° data differ 
by 20%: 364 ±7 GeV -2 for the exclusive p versus 
303 ± 10 GeV~ 2 for the Coulomb breakup events. 

The different b values may be attributed to the 
different impact parameter distributions caused 
by the nuclear breakup tagging in the minimum 
bias data. The photon flux depends on the im- 
pact parameter as 1/& 2 . When b ss fewi?^, p are 
more likely to be produced on the side of the tar- 
get near the photon emitter than on the far side. 
p° production is concentrated on the near side, 
leading to a smaller effective production volume 
and the smaller b. This near-side skewing affects 
the interference slightly, but is not included in 
current calculations. 

Systematic errors come from a variety of 
sources. We have studied the effect of detector 
distortions (primarily px smearing) by turning off 
the detector simulation and comparing raw sim- 
ulations with reconstructed data. This lowered 
c by 18%; if the detector simulation is 75% cor- 
rect (a very conservative assumption), then de- 
tector effects are a less than 5% systematic un- 




Figure 4. Efficiency corrected t± spectrum for p° from (a) mutual dissociation with 0.1 < \y\ < 0.5, 
(b) mutual dissociation with 0.5 < \y\ < 1.0, (c) topology trigger with 0.1 < \y\ < 0.5 and (d) topology 
trigger with 0.5 < \y\ < 1.0. The points are the data, while the solid line is the result of the fit to Eq. 4. 
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certainty. Minor systematic uncertainties come 
from backgrounds and the fitting procedure. We 
have compared the data and simulations for a va- 
riety of kinematic and detector-based variables, 
and found good agreement. We estimate a pre- 
liminary experimental systematic uncertainty of 
8%. 

This analysis depends on the shape of R(t), 
which is calculated following Ref. [8]. For the 
broader picture, it is necessary to consider the 
simplifications in Ref. [8]. The calculation aver- 
ages the photon flux over the surface of the nu- 
cleus, rather than incorporating the proper 1/b 2 
weighting. This may partly explain the relatively 
poor fit for the large- \y\ Coulomb breakup fit. 
We estimate that the effective impact parameter 
should be within 10% of the actual ion-ion separa- 
tion for the coulomb breakup data, and 3% for the 
exclusive p. The uncertainties in the calculations 
should be at most a 15% effect. With this, the in- 
terference is 93±6(stat.) ±8(syst.) ±15(theory)% 
(STAR preliminary) of that expected. 

The p decays rapidly, with "fficr <C (b) and the 
two p decay points are well separated in space- 
time. In the usual space-time picture, the decays 
occur independently, and any interference must 
involve the final state ir + ir~. Interference must 
involve identical final states from the two sources. 
However, given the large available phase space for 
the decays, this is very unlikely for independent 
decays. 

One possible interpretation of this result is 
given in Ref. [10]. In it, the interference occurs 
because the post-decay wave function includes 
amplitudes for all possible final states, then the 
amplitudes for identical states subtract, and the 
interference is visible. Because of the two sources, 
the 7r + 7r~ wave function is non-factorizable, and 
thus exhibits the Einstein-Podolsky-Rosen para- 
dox [25]. We find that the decoherence, 1 — c, due 
to environmental or other factors is less than 43% 
at the 90% confidence level. 

5. e + e~ Pairs 

The cross section for the production of e+e~ 
pairs is peaked near threshold, with pair mass 
M ee « 3.5m, where m is the electron mass. The 



pairs are produced predominantly in a forward- 
backward topology, with the produced with 
large longitudinal (along the beampipc) momen- 
tum and small pr- This complicates experimen- 
tal detection at an ion collider, and reduces the 
statistics available to any analysis. This analysis 
used data about 800,000 minimum bias triggers 
taken in a 0.25 T magnetic field. The lower mag- 
netic field (compared to the 0.5 T full field data) 
increased the sensitivity to low px particles. The 
minimum bias trigger freed the electrons from the 
requirement that they have a high enough px to 
reach the CTB. 

This analysis [26] used tracks with px > 65 
MeV/c and pseudorapidity \t]\ < 1.15. In this re- 
gion, tracking efficiency for electrons was above 
80%. The electrons were identified by dE/dx. To 
get good electron/hadron separation, we required 
the tracks to have p < 130 MeV/c. In this mo- 
mentum region, electrons had considerably higher 
dE/dx than hadrons, and the identification effi- 
ciency was nearly 100%, with minimal contamina- 
tion. The cross section falls steeply with increas- 
ing M ee , so few leptons were expected with higher 
momenta. We also required the events to have 
pair pt < 100 MeV/c, pair rapidity \y\ < 1.15 
and pair mass 140 MeV < M ee < 265 McV. These 
cuts selected a sample of 52 events. 

The major backgrounds in this analysis are due 
to misidentified coherent 7r + 7r~ pairs and inco- 
herent hadronic backgrounds. The former can 
be calculated from the known cross section and 
misidentification probability. The latter can be 
studied by using like sign (tt + tt + ) pairs. The total 
background is estimated to be about 1.1 events. 

The integrated luminosity of 94 ± 9 mb _1 was 
estimated by counting hadronic events with at 
least 8 tracks. This selected 80% of all hadronic 
gold-gold interactions. The integrated luminosity 
is based on an assumed gold-gold hadronic cross 
section of 7.2 barns. 

Figure 5 shows the cross section as functions 
of pair mass and pr- The data is compared with 
two calculations. The first, shown by the solid 
line, uses Eq. 3, with a equivalent photon calcu- 
lation for P ee (b). The second, provided by Kai 
Hencken, was a quantum electrodynamic calcula- 
tion, including photon virtuality. The major dif- 
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Figure 5. Differential Cross sections for pair (a) mass and (b) pr- The data (points) are compared with 
the results of a equivalent photon calculation (solid) and a quantum electrodynamics calculation (dashed 
line). All data and calculations are restricted to the kinematic region described in the text. 



ference between the results of the two calculations 
occurs for pr < 15 MeV/c. The data strongly fa- 
vors the QED calculation, by 4a. 

As expected, the {pr) for e + e~ production is 
smaller than for p° production. This is because 
the Pomeron pj- scale is roughly H/Ra, while the 
photon pt scale is very roughly M ee /j w Ti/b. 

Within the kinematic acceptance region, we 
find a cross section a = 1.6 ±0.2 ±0.3 mb (STAR 
preliminary). The major systematic errors are 
due to tracking and vertexing uncertainty, and 
the integrated luminosity. This is 1.2a lower than 
the equivalent photon prediction of 2.1 mb, and 
0.8a below the QED calculation of aQED = 1.9 
mb for the same kinematic cuts. This measure- 
ment and the QED calculation can be used to put 
limits on possible higher order corrections Aa to 
the cross section. At a 90% confidence level, we 
find -0.5a QE D < Act < 0.2a QED . 

6. Conclusions 

STAR has studied a number of photonuclear 
and two-photon reactions. We observe both ex- 
clusive p° production and p° production accom- 
panied by nuclear excitation. The cross sec- 
tions and rapidity distributions match the pre- 
dictions of the soft Pomeron model. We observe 



destructive interference between the two produc- 
tion sites, at the expected level, and set an upper 
limit on decoherence. Finally, we observe e + e~ 
pair production accompanied by mutual Coulomb 
dissociation. The kinematic distributions match 
those predicted by lowest order quantum electro- 
dynamics. 

We thank Tony Baltz for providing the calcu- 
lations of nuclear breakup and Kai Hencken for 
providing the QED calculations. This work was 
supported by the U.S. DOE under contract num- 
ber DE-AC-03-76SF00098. 
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